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Abstract 

In rhombohedral CaMn70i2, an improper ferroelectric polarization of magnitude 2870 /iC m~^ 
is induced by an incommensurate helical magnetic structure that evolves below Tni = 90 K. The 
electric polarization was found to be constrained to the high symmetry three-fold rotation axis 
of the crystal structure, perpendicular to the in-plane rotation of the magnetic moments. The 
multiferroicity is explained by the ferroaxial coupling mechanism, which in CaMn70i2 gives rise 
to the largest magnetically induced, electric polarization measured to date. 

PACS numbers: 75.85. +t, 75.25.-j, 77.80.-e 
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The coexistence of electric polarization and long-range magnetic order in single phase, 
multiferroic materials has been the topic of considerable research. The ultimate goal is to find 
a compound that displays a large improper ferroelectric polarization strongly coupled to the 
magnetic structure at high temperatures, and hence suitable for application in technology. A 
direct coupling between magnetism and ferroelectricity has been demonstrated in a number 
of multiferroics, for example TbMnOg fl, TbMnsOg Q, NigVaOg 3], and MnW04 E], 
however the electric polarization tends to be small (~ 10^ fiC m~^) and the Neel temperature 
ow (< 40 K). Few magneto-electric multiferroics, for example the hexaferrites sl and CuO 
6|, are polar close to room temperature (RT), and no such material has shown a large 
polarization comparable to that of proper ferroelectrics. 

Here, we determine the origin of the magnetically induced ferroelectricity in CaMn70i2 
0] through the analysis of electric polarization, magnetic susceptibility, and neutron powder 
diffraction measurements. We show that the magnetically induced polarization, the largest 
measured to date in any compound, develops at 90 K, oriented parallel to the c axis. The 
improper ferroelectricity is induced by a helical magnetic structure with incommensurate 
propagation vector fci = (0, 1, 0.963). Contrary to conventional magnetic multiferroics, the 
electric polarization was found to lie perpendicular to the spin rotation plane, and is allowed 
through coupling to the ferroaxial component of the crystal structure |8|. 

High quality single crystals of CaMn70i2 with approximate dimensions 300 x 300 x 300 /im 
were grown by the flux method. A 5 g mixture of CaCl2:Mn02 (weight ratio 1:3) was heated 
for 24 h in an alumina crucible in air at 850 °C, and then cooled at a rate of 5° h^^. 
X-ray diffraction and bulk property measurements showed that the single crystals were 
not twinned at RT. A 1 g powder sample was prepared by grinding single crystals and 
sieving through a 35 fim mesh. The electric polarization was determined by integrating 
the pyroelectric current measured on warming at a rate of 1 K min^^, having field cooled 
the sample in 4.4 kV cm~^. Magnetic susceptibility measurements were performed using a 
Quantum Design, Magnetic Properties Measurement System. Neutron powder diffraction 
(NPD) data were collected on the WISH time-of-fiight diffractometer at ISIS, UK j^. Data 
were collected at RT, and on warming between 1.5 K and 98 K in 3 K steps using a helium 
cryostat. The sample was enclosed in a cylindrical vanadium can. Magnetic structure 



refinements were performed using FullProf 10|] against data measured in detector banks at 



average 29 values of 58°, 90°, 122°, and 154°, each covering 32° of the scattering plane. 




FIG. 1. (Color online) The crystal structure of CaMn70i2 in both hexagonal and cubic lattices. 
Calcium ions are shown in yellow, Mnl in red, Mn2 in blue, Mn3 in green, and oxygen in black. 



At high temperatures, CaMn70i2 crystalhzes into a distorted perovskite structure with 
cubic space group ImS [ll| (see Fig. [1]). The crystal lattice undergoes a rhombohedral 
distortion on cooling through a first order phase transition at ~440 K, giving the RT space 



group R3 [12|, M with lattice vectors (—1, 1, 0), (0, —1, 1), and (|, |, |) with respect 
to those of the cubic unit cell. Here, all real and reciprocal space vectors are given in the 
hexagonal basis. In R3 there are three symmetry inequivalent manganese sites that we label 
Mnl, Mn2, and Mn3. The Mnl ions occupy sites of symmetry 1 (Wyckoff position 9e) and 
form a 3:1 ordered occupation of the pseudo-cubic perovskite A-sites (the minority A-sites 
are occupied by non-magnetic calcium ions). Mn2 and Mn3 ions occupy sites of symmetry 1 
and 3 (Wyckoff positions 9d and 36), respectively. They are octahedrally coordinated with 
oxygen and correspond to the pseudo-cubic B-sites with a ratio of three Mn2 to one Mn3. 
At RT Mnl ions have a valence of +3, and the Mn2 and Mn3 ions are mixed- valent (+3.25). 
A further isostructural charge-ordering transition at 250 K leaves nominal valences of +3 
and +4 on the Mn2 and Mn3 ions, respectively 14 1. 

CaMn70i2 undergoes a further two phase transitions, Tm = 90 K and Tn2 = 48 K, both 
of which have magnetic origin [isl . IGj. Measurements of polycrystalline samples showed 
that the development of long-range antiferromagnetic order below Tni coincides with the 
onset of a large ferroelectric polarization 0]. In the present study, the electric polarization 
was measured on our single crystals perpendicular to the natural facets (parallel to the 
<100> pseudo-cubic axes) and along the three- fold axis of the rhombohedral cell, as shown 
in Fig. [2^. The three measurements along (0, 2, 1), (2, 0, -1), and (2, -2, 1) showed 
polarization of equal magnitude within systematic error. In rhombohedral symmetry this 
equivalence may arise through one of two scenarios; either the electric polarization lies 



within the ab plane, forming a 120° domain configuration, or the polarization is parallel 
to the three-fold c-axis. The fourth measurement along (0, 0, 1) conclusively demonstrates 
that the polarization vector, P, lies parallel to the three-fold axis, where the polarization 
components along <0, 2, 1> are reduced by a factor l/VS. Remarkably, the saturation 
value |P| = 2870 fiC m~^ is the largest measured for magnetically-driven ferroelectrics to 
date, and is at least four times larger than that of the prototypal multiferroic, TbMnOa 

The magnetic susceptibilities measured parallel and perpendicular to the c axis (Fig. [2]d) 
display accidents coinciding with T^i and Tn2- At T^i, a small downturn of the in-plane 
susceptibility suggests antiferromagnetic ordering of the Mn moments in the hexagonal basal 
plane. At T^2, the behavior is more complex with a small peak in the c axis susceptibility, 
possibly indicative of an additional ordering along this direction, together with a rise of the 
in-plane susceptibility. 

Further insights into the nature of the magnetic state of CaMn70i2 were gained through 
NPD measurements. A large number (>15) of magnetic refiections were observed be- 
low Tni (Fig. [3]) indicative of long-range magnetic ordering. In the temperature range 
< T < Tni, all magnetic refiections can be indexed with the propagation vector 
fci = (0, 1, 0.963) (phase AFM I). Below (AFM II) however, it is impossible to index 
all magnetic refiections with a single A;-vector. Instead, one observes the coexistence of two 
modulations along the same symmetry line of the Brillouin zone, with respective propagation 
vectors ^2 = (0, 1, 0.880) and ^3 = (0, 1, 1.042) at 5 K. In fact, all three propagation vec- 
tors lie along the same {hc,hc,hc) pseudo-cubic line of symmetry. The thermodiffractogram 
clearly shows that the magnetic peaks indexed by ki above T^i split symmetrically at T^2 
into the two contributions corresponding to fc2=fci-(0,0,5) and k3=ki+{0,0,6), as shown in 
Fig. 121:;, as well as in the temperature dependence of the integrated intensities (Fig. |2ll). We 
note that this magnetic phase diagram contradicts earlier neutron diffraction results which 
claimed a ferrimagnetic component in AFM I, incompatible with the bulk magnetization 



measurements 15 1 



In the following we focus on the magnetic structure of the AFM I phase, which is key 
to understanding the onset of the large electric polarization in CaMn70i2. The magnetic 
structure of AFM II is more complex due to the presence of two A;-vectors and will require 
a single crystal, neutron study when larger specimens are available. For AFM I, symmetry 
analysis indicates that there are three irreducible representations (irreps) in the little group. 




FIG. 2. (Color online) a) The electric polarization of CaMn70i2, measured parallel to the hexag- 
onal c-axis and the three pseudo-cubic <100> axes. Inset: A small increase in P at Tn2- b) 
The magnetic susceptibility parallel and perpendicular to the hexagonal c-axis, measured with 
H = 500 Oe, in both zero field cooled and field cooled (500 Oe) conditions. Inset: Reduction of 
the in-plane magnetic susceptibility at Tni. c) The variation of the incommensurate magnetic 
propagation vectors ki (black circles), ^2 (blue triangles), and fcs (red squares), refined against 
neutron powder diffraction data. The error bars are within the size of the data point, d) The 
integrated intensity of the (1, -1, 1-kz) magnetic reflection at d ~ oA. The black circles, blue 
triangles and red squares correspond to fci, k2, and fca, respectively. 
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FIG. 3. (Color online) a) Rietveld refinement against neutron powder diffraction data at 65 K. b) 
Le Bail fit to data measured at 5 K. Top and bottom tick marks indicate the nuclear and magnetic 
peaks, respectively. The measured and calculated profiles are shown with dots and a continuous 
line, respectively. A difference curve (observed— calculated) is shown at the bottom. 

all one-dimensional. All irreps appear in the decomposition of the magnetic representation 
for each site, suggesting that all Mn moments can be simultaneously ordered even if a single 
irrep is involved at the phase transition. For Mnl and Mn2 sites, the possible magnetic 
arrangements spanned by the three irreps correspond to the symmetry-adapted modes of 
a system of spins on a triangle, as displayed in Fig. Ht. For a structure with equal mo- 
ments, the Fi and F2 solutions lead to the extinction of the (0, 0, l+kz) reflection due to 
the vanishing net moment in the triangle (null magnetic structure factor). However, such 
peaks are observed in the AFM I diffraction data (labelled in Fig. dictating that the 
magnetic structure has F3 symmetry, corresponding to a ferromagnetic arrangement of both 
Mnl and Mn2 triangles, as shown in Fig. Ht. Refinement of the data at 65 K indicated 
that the moments are confined to lie in the ab plane. The Mn^+ moments follow a helical 
modulation along the c axis with a circular envelop, imposed by symmetry. The Mn'^"^ also 
follow a circular, helical modulation along c, since a refinement with an elliptical envelop 
(allowed by symmetry) did not lead to any improvement of the fit. The final refinement 
was conducted with five free parameters; the moment magnitudes of the three inequivalent 
manganese sites and the relative phases between each helix. The refinement is shown in Fig. 
|3^ and the results are given in Table [11 Reliability factors of Rp = 6.65% and R^p = 5.62% 
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FIG. 4. (Color online) (a), (b) The magnetic structure of CaMn70i2 in the ab and ac planes, 
respectively. Mnl, Mn2, and Mn3 ions are shown in red, black and yellow, respectively. The 
moments rotate in the a6-plane with a circular envelope as depicted in (b), perpendicular to the 
electric polarization || c-axis. c) The three triangular modes of symmetry equivalent Mnl or Mn2 
ions that correspond to Fi, F2, and F3. The magnetic structure of CaMn70i2 has F3 symmetry, 
and the respective triangular mode of equivalent Mnl sites is shown by the heavy line in (a). 

were achieved for the 58° data, with a global (all banks) = 12.0. 

The global magnetic structure can be described as follows: Mn^"*" ions form a buckled 
kagome lattice, as shown in Fig. H^, in which the spin arrangement follows the so-called 
(■\/3,\/3)^ structure, one of the two possible ordered states of the classical kagome antiferro- 
magnet [17]. Each triangle has a 120° spin configuration. The buckled kagome layers are 
then stacked along the c axis (Fig. Hb), with a long period modulation (A;2=0.963) corre- 
sponding to a small global rotation of ~ 4°. Mn^"*" ions occupy the hexagonal voids of the 
kagome lattice and their moments are aligned ferromagnetically with respect to their nearest 
neighbor Mn^"*" ions in each layer. This chiral magnetic structure breaks inversion symmetry 
with magnetic point group 31'. 
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TABLE I. Magnetic structure parameters at 65 K 



Atom Valence x y z (/iR) 4> (tt rad) 

Mnl +3 0.5 2.31(3) 

Mn2 +3 0.5 0.5 2.33(3) -0.03(5) 

Mn3 +4 0.5 2.01(8) 0.84(4) 



The electric polarization was found to lie perpendicular to the spin rotation plane of the 
lelical magnetic structure. This rules out the conventional inverse Dyzaloshinskii-Moriya 
18| or spin current [19i] models for induced ferroelectricity, in which the electric polarization 
is constrained to the plane of rotation of the spins. In contrast, the magneto-electric prop- 
erties of CaMn70i2 are allowed through coupling to the ferroaxial component of the crystal 
structure Jsjl , which also explains the multiferroic properties of other helical magnetic struc- 



tures 
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-22l|. The nonpolar paramagnetic phase belongs to the 3 ferroaxial crystal class, for 
which a homogeneous structural rotation exists with respect to the high-temperature cubic 
group, represented by an axial vector A. Symmetry considerations show that the chirality 
of the helical magnetic structure, a = ■ {SixSj), a time-even, parity-odd pseudoscalar, 
couples to A and the polarization vector (necessarily || A), P^, to form the invariant '-fPzCrA, 
where 7 is a coupling constant. Considering the dielectric term in the free energy expan- 
sion, —Pli the stability condition leads to P^oc -e^^crA. In the first approximation, A is 
constant at the temperatures in question (well below the critical temperature of the struc- 
tural transition), therefore P z varies linearly with a. In the mean- field approximation the 
magnetic critical exponent is 0.5. a, and hence P z-, should vary as the square of the magne- 
tization, {SixSj), and therefore linearly with temperature, as observed in the temperature 
dependence of the electrical polarization in CaMn70i2 (except very close to the transition). 

The exact microscopic coupling between spins and electric dipoles is unclear due to the 
complexity associated with the numerous competing exchange interactions. However, it is 
likely that exchange-striction coupled to the Mn^"*" / Mn^"*" charge order is responsible for 
such a large electric polarization. 

To conclude, we have found that in CaMn70i2 a giant improper ferroelectric polarization 
is induced parallel to the c axis by an in-plane helical structure at Tni = 90 K. The 
multiferroicity can be explained by the ferroaxial coupling of the magnetic chirality to the 
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macroscopic structural rotations associated with the -R3 paramagnetic phase. 

The work done at the University of Oxford was funded by an EPSRC grant, number 
EP/J003557/1, entitled "New Concepts in Multiferroics and Magnetoelectrics" . 
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